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Abstract. The zero locus of a function / on a graph G is de¬ 
fined as the graph for which the vertex set consists of all complete 
subgraphs of G, on which / changes sign and where x, y are con¬ 
nected if one is contained in the other. For d-graphs, finite simple 
graphs for which every unit sphere is a d-sphere, the zero locus of 
(/ — c) is a (d — l)-graph for all c different from the range of /. If 
this Sard lemma is inductively applied to an ordered list functions 
/i,..., /fc in which the functions are extended on the level surfaces, 
the set of critical values (ci,..., c^) for which E — c = 0 is not a 
(d — fc)-graph is a finite set. This discrete Sard result allows to 
construct explicit graphs triangulating a given algebraic set. We 
also look at a second setup: for a function F from the vertex set 
to R^, we give conditions for which the simultaneous discrete al¬ 
gebraic set {F = c} defined as the set of simplices of dimension 
€ {k,k + 1,... ,n} on which all ft change sign, is a (d — fc)-graph 
in the barycentric refinement of G. This maximal rank condi¬ 
tion is adapted from the continuum and the graph {F = c} is a 
(n — A:)-graph. While now, the critical values can have positive 
measure, we are closer to calculus: for k = 2 for example, ex¬ 
trema of functions / under a constraint {g = c} happen at points, 
where the gradients of / and g are parallel V/ = AV5, the La¬ 
grange equations on the discrete network. As for an application, 
we illustrate eigenfunctions of d-graphs and especially the second 
eigenvector of 3-spheres, which by Courant-Fiedler has exactly two 
nodal regions. The separating nodal surface of the second eigen¬ 
function fi always appears to be a 2-sphere in experiments. By 
Jordan-Schoenfliess, both nodal regions would then be 3-balls and 
the double nodal curve {/i = 0, /2 = 0} would be an un-knotted 
curve in the 3-sphere. Graph theory allows to approach such unex¬ 
plored concepts experimentally, as the corresponding question are 
open even classically for nodal surfaces of the ground state of the 
Laplacian of a Riemannian 3-sphere M. 
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1. Introduction 

We explore vector-valued functions F on the vertex set of a finite sim¬ 
ple graph G = (V, E). Most of the notions introduced here are defined 
for general hnite simple graphs. But as we are interested in Lagrange 
extremization, Morse and Sard type results in graph theory as well 
as questions in the spectral theory of the Laplacian on graphs related 
to Laplacians of Riemannian manifolds, we often assume G to be a 
d-graph, which is a finite simple graph, for which all unit spheres are 
(d — l)-spheres in the sense of Evako [20]. In a first setup, more suited 
for Sard, for all except finitely many choices of c G the graph 
{F = c} = {/i = Cl,..., fk = Ck} is a (d — fc)-graph, in a second setup, 
closer to classical calculus, we need to satisfy locally a maximal rank 
condition to assure that the graph representing the discrete algebraic 
set {F = c} is a (d — fc)-graph. 

The first part of the story parallels classical calculus and deals with the 
concept of level surface. It starts with a pleasant surprise when looking 
at a single function: there is a strong Sard regularity for a d-graph G: 
given a hyper-surface Gf{x) = {/ = c} defined as the graph with ver¬ 
tex set consisting of all simplices on which / — c changes sign, the graph 
{/ = c} is a (d — l)-graph if c is not a value taken by /. The topology 
of Gf{c) changes only for parameter values c contained in the hnite set 
f{V). This observation was obtained when studying coloring problems 
HailHlDE], as a locally injective function on the vertex set is the same 
than a vertex coloring. Before, in P E], we just looked at the edges, on 
which / changes sign and then completed the graph artihcially. Now, 
we have this automatic. We hope to apply this to investigate nodal 
regions of eigenfunctions of the Laplacian of a graph, where we believe 
the answers to be the same for compact d-dimensional Riemannian 
manifolds or hnite d-graphs. With the context of level surfaces one has 
the opportunity to look at nodal surfaces of eigenfunctions, which are 
also known as Chladni surfaces bounding nodal regions. 

With more than one function, the situation changes as the singular set 
typically becomes larger in the discrete. This is where the story splits. 
In the commutative setup, where we look at the zero locus of all func¬ 
tions simultaneously, Sard fails, while in the setup, where an ordered 
set of functions is considered, Sard will be true: as in the classical Sard 
theorem [221123] . the set of critical values has zero measure. The differ¬ 
ence already is apparent if we take two random functions on a discrete 
3-sphere for example. Then simultaneous level set {/ = 0,5^ = 0 } is 
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a graph without triangles but it is rarely a 1 -graph, a hnite union of 
circular graphs. The reason is that the tangent space is a hnite set 
and the probability of having two parallel gradients at a vertex does 
not have zero probability. However, if we look at 5 ^ = 0 on the two 
dimensional surface / = 0, we get a hnite union of cyclic graphs. In 
general, we salvage regularity and Sard by dehning the algebraic set 
{F = c} in a diherent way by recursively building hyper surfaces: start 
with the hyper surface {/i = Ci}, then extend /2 to the new graph 
and look at {/2 = C 2 } inside {/i = Ci}. Now a vector value c E is 
always regular if c is not in the image of F applied to recursively de- 
hned graph. This is the discrete analog of the multi-dimensional Sard 
theorem in classical analysis. The order in which the functions are 
taken, matters in the discrete. But this is not a surprise as we rehne in 
each step the graph and therefore have to extend the functions to the 
barycentric rehnements. Its only for sufficiently smooth functions like 
eigenfunctions of low energy eigenvalues of the graph that the answer 
can be expected to be independent of the ordering. 

The graph theoretical approach is useful for making experiments: take 
a 3-graph G = (V, E) for example and take two real valued functions 
/, g on the vertex set V. If c is not in the image of /, we can look at 
the level surface / = c and extend 5 ^ to a function there (vertices are 
now simplices and we just average the function value g to extend the 
function). If d is not in the image of g, then the 1-graph Z = {g = d} is 
a subgraph of {/ = c}. It is a hnite set of closed curves in G 2 . In other 
words, each connected component is a knot in the 3-sphere. Unlike in 
the continuum, we do not have to worry about cases, where the knot 
intersects or self intersects. The Sard theorem assures that this will 
never happen in the discrete. We only have to assure that the value c 
is not in fiV) and d is not in giV). 

The second setup, which dehnes discrete algebraic sets Z = {fi = 
Cl,fk = Ck} in the Barycentric rehnement Gi of G requires a few 
dehnitions. We need conditions under which these sets are nonsingular 
in the sense that they again form a (d — A;)-graph, graphs for which the 
unit sphere S{x) is a sphere at every vertex. In general, the situation 
is the same as in the continuum, where varieties are not necessarily 
manifolds. The dehnition is straightforward: dehne {F = c } as the 
set of simplices of dimension in {k, ..., n} in G on which all functions 
fj change sign simultaneously. The graph {F = c} is a subgraph of 
the Barycentric rehnement Gi = G x 1 of G. Conditions for regu¬ 
larity could be formulated locally in terms of spheres in spheres. In 
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Figure 1. We see an example of a level surface 
^/(c) = {/ = ^ 3-sphere G. The graph G is a 

Barycentric rehnement of a 3D octahedron with 8 ver¬ 
tices with simplex cardinalities v = (80,464,768,384) 
satisfying x{G) = 1)*D = 0. The zero Euler char¬ 

acteristic is a property shared by all 3-graphs like 3- 
manifolds. The graph G admits an involution T such 
that G/T is a discrete projective 3-space. The alge¬ 
braic level set G/(c) displayed here is a 2-graph, a two- 
dimensional surface within the Barycentric refinement 
Gi of G. The ambient sphere Gi already has 1696 ver¬ 
tices and 10912 edges and 18432 triangles, 9216 tetra- 
hedra. The surface {/ = c} is a discrete analogue of 
a 2-dimensional projective variety. By the Sard lemma, 
the level sets Gf{c) = {/ = c} are always nonsingular if 
c ^ f{V). Furthermore the surface Gf{c) is 4-colorable. 
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some sense, a locally projective tangent bnndle with discrete projec¬ 
tive spaces associated to nnit spheres replaces the tangent bnndle. It 
is still possible to dehne a vector space structnre at each point and 
dehne a discrete gradient V/ of a fnnction. If G is a d-graph, x is 
a complete Kd+i snbgraph, then Vf{x,xo) is dehned as the vector 
V/ = (sign(/(a;i) - /(xq)), ... ,sign(/(xd) - /(xq))). The local injec¬ 
tivity condition means that V/(x,Xo) is not zero for all d-simplices 
containing xq. This leads to the Sard Lemma Q. Bnt now, when 
changing c, we also need to look at the valnes c, = /(xj) taken by the 
fnnction / and investigate whether they are critical. What happens if c 
passes a valne q = /(xj)? If / < q is homotopic to / < q, then the set 
{/ = c} n S{x) is a (d — l)-sphere which by Jordan-Brower-Schoenflies 
[20] divides S'(x) into two complementary balls ^“(x), S~^{x). The 
Poincare-Hopf index 1 — x(5'“(x)) = f/(x) [7] is then zero. At a local 
minimnm of / for example, the graph S~ (x) is empty so that the index 
is 1. At a local maximnm, S~^{x) is empty and the index is either 1 or 
— 1 depending on whether the dimension of G is even or odd. In two 
dimensions, where we look at discrete mnlti-variable calcnlus, maxima 
and minima have index 1 and saddle points have index —1. 

The set of hypersnrfaces Mj = {/ = cj} form a contour map for 
which the individual leaves in general have different homotopy types. 
Topological transitions can happen only at values f{V) of / on the 
vertex set V of the graph. We can relate the symmetric index jf{x) = 
[if{x) -fi_/(x )]/2 at a vertex x of G with the (d — 2 )-dimensional graph 
{/ = Cfc} in the unit sphere S'(x). We have called this the graph Bf{x) 
and showed that the graph can be completed to become a d-graph. 
This completion can be done now more elegantly. As pointed out in 
[13] . for d = 4 for example, we get for every locally injective function 
/ a 2-dimensional surface i?/, the disjoint union of Bf{x). We see that 
a function / not only dehnes (d — l)-graphs G/(c) = {/ = c} but also 
(d — 2 )-graphs i?/(x) called “central surfaces” for every vertex x. 

If we have more than one function as constraints, the singularity struc¬ 
ture of {/i = Cl,..., /fc = Cfc } is more complicated and resembles the 
classical situation, where singularities can occur. Also the Lagrange 
setup, where we mazimize or minimize functions under constraints, is 
very similar to the classical situation. The higher complexity entering 
with 2 or more functions is no surprise as also for classical algebraic sets 
dehned as the zero locus of hnitely many polynomials, the case of sev¬ 
eral functions is harder to analyze. In calculus, when studying extrema 


6 


OLIVER KNILL 


of a function / under constraints g, following Lagrange, one is inter¬ 
ested in critical points of / and g as well as places, where the gradients 
are parallel. Lets assume that we have two functions /, g on the vertex 
sets of a geometric graph. The intersection {f = c,g = d] with S{x) 
is then a sphere of co-dimension 2. If we have 2 constraints and G is 
4-dimensional for example, then S{x) is a 3-sphere and S{x) fl {F = 0} 
is a knot inside S{x). This can already be complicated. By triangulat¬ 
ing Seifert surfaces associated by a knot, one can see that any knot can 
occur as a co-dimension 2 curve of a graph. Briskorn sphere examples 
allow to make this explicit in the case of torus knots. With more than 
one function, Poincare-Hopf indices form a discrete 1-form valued grid 
because changing any of the cj can change the Euler characteristic. 
This allows to express the Euler characteristic as a discrete line inte¬ 
gral in the discrete target set of F. There are now many Poincare-Hopf 
theorems, for every deformation path, there is one. 

If we look at the set F = c, where all functions change sign, singular 
values for F = (/i,..., fk) are vectors c = (ci ,... ,Ck) in for which 
the graph {F = c} is not a (d — /c)-graph or values c = F{x) taken on 
by F on vertices x. Lets look at the very special example, where all 
fj are the same function. Now, near the diagonal q = c, there is an 
entire neighborhood of parameter values, where regularity fails. We see 
that the Sard statement does not hold in this commutative setting. We 
therefore also look at the non-commutative setup, where we £x /i = ci 
first, then look at the level surface {/2 = C 2 } on the surface {/i = Ci} 
and proceed inductively. Sard is now more obvious, but the sets F = 0 
depend on the order, in which the functions fj have been chosen. This 
order dependence is no surprise in a quantum setting, if we look at fj 
as observables. It simply depends in which order we measure and fix 
the fj. 

There are analogies to Morse theory in the continuum; in the case of 
one single function, we can single out a nice class of graphs and func¬ 
tions, which lead to a discrete Morse theory. The geometric graph G as 
well as functions on vertices are the only ingredients. One can now as¬ 
sign a Morse index m{x) at a critical point and have if{x) = (— 

The requirement is an adaptation of a reformulation of the definition of 
being Morse means in the continuum that for a small sphere S{x), the 
set S{x) n {f{y) = f{x)} is a product Sm-i x Sd-i-m- For example for 
m = l,d = 2, we have a saddle point, where S{x) intersects the level 
surface in 4 points forming Sq x S'q. In the enhanced picture, where 
we look at the graph Gi, we have this regularity more likely. One can 
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extend the function to the new graph and repeat until one get a Morse 
function. 

We can use graphs described by hnitely many equations in order to 
construct examples of graphs to illustrate classical calculus like Stokes 
or Gauss theorem or surfaces to illustrate classical multivariable cal¬ 
culus. The notion of d-graphs has evolved from P El 1131 HI] to [IH], 
where it reached its hnal form. While hnishing up [2D] a literature 
search showed that spheres had been dehned in a similar way already 
by Evako earlier on. 

The enhanced Barycentric rehnement graph Gi = G x Ki is a. regu¬ 
larizes graph which helps to study Jordan-Brouwer questions in graph 
theory [20| and introduce a product structure on graphs which is com¬ 
patible with cohomology [2T|. It also illustrates the Brouwer hxed 
point theorem m, as the hxed simplices on G can be seen as hxed 
points on Gi. The simplex picture is also useful as the Dirac operator 
D = d + d* on a graph builds on it imiiD]. The graph spectra of suc¬ 
cessive Barycentric rehnements converges universally, only depending 
on the size of the largest complete subgraph na. 

An other application of the present Sard analysis is a simplihcation 
of [8] which assures that the curvature K is identically zero for d = 
(2m -f l)-dimensional geometric graphs: write the symmetric index 
jf{x) = if{x) + i-f{x) in terms of the Euler characteristic of the 
{d — 2)-dimensional graph Gf{x) obtained by taking the hypersurface 
{y I f{y) = f{^)} ill the unit sphere S{x). For odd-dimensional geo¬ 
metric d-graphs, this symmetric index is is —x{Bf{x))/2 = jf{x), while 
for even-dimensional graphs, it is 1 — x{Bf{x))/2 = jf{x). One of the 
corollaries given here is that if / is locally injective, then Bf{x) is al¬ 
ways a geometric (d — 2)-graph if G is a d-graph. In [9], we called 
the uncompleted graphs polytopes and showed that one can complete 
them. Now, since the expectation of j/(x) is curvature, [HI ED], this 
gives an immediate proof that odd-dimensional d-graphs have constant 
0 curvature. Zero curvature follows for odd-dimensional graphs. The 
regularity allows to interpret Euler curvature as an average of two- 
dimensional sectional curvatures. Euler characteristic written as the 
expectation of these sectional curvature averages is close to Hilbert ac¬ 
tion. It shows that the quantized functional “Euler characteristic” is 
not only geometrically relevant but that it has physical potential. 


OLIVER KNILL 


2. Level surfaces 

The study of level surfaces {/ = c} in a graph is not only part of discrete 
differential topology. It belongs already to discretized multivariable 
calculus, where surfaces {f{x,y,z) = 0} in space or curves {f{x,y) = 
0} in the plane are central objects. Some would call calculus on graphs 
“quantum calculus”. We want to understand under which conditions a 
sequence of k real valued functions F = (/i,..., /fc) on the vertex set 
of a graph leads to co-dimension k graphs {F = 0}. We start with the 
case k = 1, where we have a level surface {/ = c}. 

Definition 1. Given a finite simple graph G = (V, E), define the level 
hyper surface {/ = c} as the graph G' = (W, E') for which the vertex 
set V is the set of simplices x in G, on which the function / — c changes 
sign in the sense that there are vertices in x for which / — c < 0 and 
vertices in x for which f — c > 0. A pair of simplices (x, y) &V' x V is 
in the edge set E' of G' if x is a subgraph of ?/ or if ?/ is a subgraph of 
X. In the case c = 0, the graph {/ = 0} is also called the zero locus 
of/. 

Remarks. 

1) Instead of taking a real-valued function, we could take a function 
taking values in an ordered field. We need an ordering as we need to 
tell, where a function “changes sign”. 

2 ) Most of the time we will assume c is not a value taken by /. An 
alternative would be to include all simplices which contain a vertex on 
which / = 0. 



Figure 2. A level curve in a planar graph. A level 
surface in a 3-sphere. 
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Examples. 

1) Let G be a 2-sphere like for example an icosahedron. Let / be 1 on 
exactly one vertex x (a discrete Dirac delta fnnction) and 0 everywhere 
else. Now {/ = 1/2} consists of all edges and triangles containing x. 
They form a circnlar graph and {/ = 0} = G 2 deg(a:)- 

2) Let G = Cn x Cn X Cn he a. discrete 3-torus and let / be a function 
which is 1 on a circular closed graph and 0 else. Then {/ = 1/2} is a 
2 -dimensional torus. 


3. The Sard lemma 

The following dehnitions are recursive and were hrst put forward by 
Evako. See [20] for our hnal version. 

Definition 2. A d-sphere is a hnite simple graph for which unit sphere 
S{x) is a (d — l)-sphere and such that removing a single vertex from 
the graph renders the graph contractible. Inductively, a graph is con¬ 
tractible, if there exists a vertex such that both S{x) and the graph 
generated with vertices without x are contractible. 

Definition 3. A d-graph G is a hnite simple graph for which every 
unit sphere S{x) is a (d — l)-sphere. 

Examples. 

1 ) A 1 -graph is a hnite union of circular graphs, for which each con¬ 
nectivity component has 4 or more vertices. 

2) The icosahedron and octahedron graphs are both 2-graphs. In the 
hrst case, the unit spheres are G5, in the second case, the unit spheres 
are G4. 

3) In [19] we have classihed all Platonic d-graphs using Gauss Bon¬ 
net [5|. Inductively, a d-graph G is called Platonic, if there exists a 
(d — l)-graph H which is Platonic such that all unit spheres of G are 
isomorphic to H. In dimension d = 3, there are only two Platonic 
graphs, the 16 and 600 cell. In dimensions d > 4, only the cross poly¬ 
topes are platonic. 

The following Sard lemma shows that we do not have to check for the 
geometric condition if we look at level surfaces: it is guaranteed, as 
long as we avoid function values in f{V). Not even local injectivity is 
needed: 

Lemma 1 (Sard lemma for real valued functions). Given a function 
f :V on the vertex set of a d-graph G = (y,E). For every 

c ^ f(y), the level surface {/ = c} is either the empty graph or a 
(d — l)-graph. 
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Proof. We have to distinguish various cases, depending on the dimen¬ 
sion of X. If X is an edge, where / changes sign, then / changes sign on 
each simplex containing x. The set of these simplices is a unit sphere 
in the Barycentric rehnement Gi of G and therefore a (d — 1) sphere. 
If a; is a triangle, then there are exactly two edges contained in x, on 
which / changes sign. The sphere S{x) in {/ = c} is a suspension 
of a (d — 2)-sphere: this is the join of Sq with Sd- 2 - In general, if x 
is a complete subgraph then the unit sphere S{x) is a join of a 
{k — 2)-dimensional sphere and a (d — /c)-dimensional sphere, which is 
a (d — 2)-dimensional sphere. As each unit sphere S{x) in {/ = c} is a 
(d — 2)-sphere, the level surface {/ = c} is a (d — l)-graph. □ 

Examples. 

1) If d = 2, and / changes sign on a triangle K^, then it changes 
sign on exactly two of its edges. If / changes sign on an edge, then it 
changes sign on exactly two of its adjacent triangles. We see that the 
level surface {/ = c} is a graph for which every vertex has exactly two 
neighbors. In other words, each unit sphere is the 0-sphere. 

2) If d = 3, and / changes sign on a tetrahedron x = ^^ 4 , then there 
are two possibilities. Either / changes sign on three edges connected 
to a vertex in which case we have 3 edges and 3 triangles in the unit 
sphere of K 4 with 4 vertices. A second possibility is that / changes 
on 4 edges and 2 triangles in which case the unit sphere consists of 6 
vertices. Now look at a triangle x = K^. It is contained in exactly 
two tetrahedra and contains two edges. The unit sphere is C 4 . Finally, 
if a; = K 2 is an edge, then all triangles and tetrahedra attached to x 
form a cyclic graph of degree G 2 n where n is the number of tetrahedra 
hinging on x. 

The Sard lemma can be used for minimal colorings for which the num¬ 
ber of colorings is exactly known: 

Corollary 2. If G is a d-graph and c is not in the range of f, then 
the surface H = {f = c} is a {d — l)-graph which is d-colorable. The 
chromatic polynomial p{x) of these graphs satisfies p{d) = d\. 

Proof. To every vertex of Gi, we can attach a “dimension” which is 
the dimension of the simplex in G it came from. This dimension is the 
coloring. It remains a coloring when looking at subgraphs. □ 

Examples. 

1) For a level surface / = 0 on a d-dimensional graph, we get a (d — 1) 
graph which is (d -|- l)-colorable. For example, for d = 2, the graph 
can be colored with 3 colors. This is minimal as any triangle already 
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needs 3 colors. It implies the graph is Eulerian: the vertex degree is 
even everywhere. 

2) If is the number of colorings with minimal color of G, we can for 
every vertex x look at the index if{x) = 1 — x{Sj{x)) where Sj{x) 
is the set of vertices y on the sphere S{x), where f{y) < f{x). Given 
a geometric graph G of dimension d, then Gi is d + 1 colorable. The 
number of colorings is (d + 1)!. We can look at the list of indices which 
are possible on each point and call this the index spectrum. The set 
of vertices x where {/ < c} changes the homotopy type are called 
critical points of /. If the index is nonzero, then we have a critical 
point because the Euler characteristic is a homotopy invariant. But 
there are also critical points with zero index, as in the continuum. 
Finding extrema of / can be done by comparing the function values of 
all vertices where if{x) is not zero or more generally, where Sj{x) is 
not contractible. At a local minimum S~ f{x) is empty. 

4. The central surface 

Besides the surface G/(c) = {/ = c}, there is for every vertex x a 
central surface of co-dimension 2 which is obtained by looking at level 
surfaces Bf{x) obtained by looking at {/ = c} inside the unit sphere 
S{x). This object was introduced in [H |13] for 4-graphs G, where 
the central surface is a 2-dimensional graph, a disjoint union of 2- 
dimensional subgraphs Bf{x) of the 3-dimensional unit spheres S{x). 

Definition 4. A real-valued function / on the vertex set of a graph 
G = {V,e) is called locally injective, if f{x) ^ f{y) for {x,y) G E. 
An other word for a locally injective function is a coloring. 

Definition 5. Given a d-graph G, a locally injective function / and 
a vertex xo, dehne the central surface Bf{xo) as the level surface 
{/ = in S{xq). It is a (d — 2) graph. The graph consists of all 

simplices in S{xq) for which / takes values smaller or larger than /(xq). 

Each of these surfaces are subgraphs of their unit sphere S'(xo). We 
have one surface for each vertex xq. On each sphere S{x,y) = S'(x) fl 
S{y), we can look at the intersection of Bf{x) and Bf{y). It consists of 
all simplices in S{x, y) where both f{z) — f{y) and f{z) — f{x) change 
sign. Sometimes they can be joined together along a circle. For exam¬ 
ple, given a 3-graph G, then the union Bj of all Bf{x) consists of all 
edges and triangles so that the max and min on each larger tetrahedron 
are attained in the edge or triangle. The following dehnition was hrst 
done in j7]: 
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Definition 6. Given a finite simple graph G and a locally injective 
real-valued function / on the vertex set V, the Poincare-Hopf index 
is defined as if (x) = 1 — x{Sj{x)), where Sj{x) is generated by {y G 
S{x) I f{y) < f{x)}. The symmetric index is defined as jf{x) = 
{if{x) +i-f{x))/2. 

The Poincare-Hopf theorem [7] tells that ~ xiG). Since 

this also holds for the function —/, we have 

= x{G) . 

xev 

The following remark made in |H] expresses jf{x) as the Euler charac¬ 
teristic of a central surface, provided the graph is geometric: 


Proposition 3. Given a d-graph G and a locally injective function f. 
If Bf[x) is the central surface, then for odd d, we have 

Jfi^) = -x{Bf{x))/2 , 

for even d, we have 

Jfi^) = ^-xiBf{x))/2 . 

Proof. Bf{x) is a (d—2)-graph in S{x), which by assumption is a (d—1)- 
graph. Since / is locally injective, the function g{y) = f{y) — f[x) does 
not take the value 0 on S{x). By the Sard lemma, the graph Bf{x) is 
a (d — 2)-graph. □ 


Corollary 4. For a d-graph G with odd dimension d, the curvature 

S2{x) 




+ 


(where Sk{x) are the number of Kk+i subgraphs of the unit sphere S{x)) 
has the property that it is constant zero for every vertex x. 


Proof. The expectation jf{x) is curvature [HI US]. Since jf{x) is iden¬ 
tically zero as the Euler characteristic of an odd dimensional (d — 2)- 
graph, also curvature is identically zero. □ 


Note that in the continuum, the Euler curvature is not even defined for 
odd dimensional graphs as the definition involves a Pfaffian |2]. Having 
the value 0 in the discrete is only natural. 


Definition 7. A function / on the vertex set is called a Morse func¬ 
tion if it is locally injective and if at every critical point, there is a 
positive integer m, such that Bf{x) = {f{y) = f{x)} within S{x) is a 
product Sm-i X Sd-i-m or the empty graph if m = 0 or m = d. The 
integer m is called Morse index of the critical point x. 
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Depending on whether m is odd or even, we have x{Bf{x)) = 4 or 0 
so that the index 1 — x{Bf{x))/2 = —1 if m is odd and 1 if m is even. 
When adding a critical point, this corresponds to add a m-dimensional 
handle. It changes the Euler characteristic by (—I)’” and changes the 
m’th cohomology by 1. 


5. Lagrange 

In this section we try to follow some of the standard calculus setup 
when extremizing functions with or without constraints. But it is done 
in a discrete setting, where space is a graph. As school calculus mostly 
deals with functions of two variables, we illustrate things primarily for 
2 -dimensional graphs, even so everything can be done in any dimen¬ 
sions. 

There are three topics related to critical points in two dimensions: 
A) extremizations without constraints, B) equilibrium points of vector 
helds and C) extremization problems with constraints which are called 
Lagrange problems. In the case A), can look at extrema of a function 
/ on the vertex set of a 2-graph, in the case B) we look at equilibrium 
points of a pair F = (/, g) of functions on the vertex set of a 2-graph, 
and hnally in the case C), we look at extrema of a function / on the 
vertex set under the constraint g = c on a 2 -graph. 

Lets first look at the “second derivative test” on graphs. Recall that a 
vertex t in a graph is a critical point of a function /, if {/ < f{x)} and 
{/ < f{^)} ci-re not homotopic. This is equivalent to the statement that 
Sj{x) = {y E S{x) I f{y) < f{x)} is a graph which is not contractible. 
The analogue of the discriminant D is the PoincarAHopf index if{x) = 
1 — x{Sj{x)). Here is the analogue of the second derivative test: 

Proposition 5 (Second derivative test). Let G he a d-graph and as¬ 
sume f is locally injective and x is a critical point. There are three 
possibilities: 

a) If Sj{x) is a {d — 1)-sphere, then x is a local maximum. 

b) If if {x) is positive and Sj{x) is empty then x is a local minimum. 

c) If d = 2 and ifif{x) is negative, then x is a type of saddle point. 

Proof. For a 2-dimensional graph, the index is nonzero if and only if 
T is a critical point because a subgraph of a circular graph has Euler 
characteristic 1 if and only if it is a contractible graph. It has Eu¬ 
ler characteristic 0 if and only if it is either empty or the full circular 
graph. In all other cases of subgraphs of a circular graph, the Euler 
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characteristic counts the number of connectivity components. 

In higher dimensions, there are cases of graphs Sj{x) having Euler 
characteristic 1 but not being contractible. In higher dimensions, Sj{x) 
can have negative Euler characteristic so that the index if{x) can be¬ 
come larger than 1. □ 

Example. 

1) For d = 2, the standard saddle point is if{x) = —1. The function 
/ changes sign on 4 points. A discrete “Monkey saddle” has index 
if{x) = —2. It is obtained for example at a vertex x for which the unit 
ball is a wheel graph with Cq boundary such that f{y) is alternating 
smaller or bigger than f{x). 

2) If d is odd and / has a local maximum, if{x) = — 1. This is analogue 
to the continuum, where D is the determinant of the Hessian. 

For simplicity, we restrict to a simple 2-dimensional situation, where we 
have two functions /, g on the vertex set V of a 2-graph G = (H, E). We 
can think of F = (/, g) it as a vector held and see the equilibrium points 
are the intersection of null-clines as in the continuum. Classically, 
the critical points of / under the constraint g = c are the places, 
where these null-clines are tangent or are degenerate in that one of the 
gradients is zero. 

Definition 8. Let F = (/, g) be two functions on the vertex set of a 
2-graph. The set of equilibrium points {{x,y) \ F{x,y) = (0,0) } is 
a zero-dimensional graph given by the set of triangles, where / and g 
both change sign. 

In more generality we have dehned the set Gp = {/i = 0,..., /*, = 
0 } as the graph whose vertex set consists of the set of simplices of 
dimension in {k,... ,d} on which all functions fj change sign and where 
two simplices are connected, if one is contained in the other. It is 
possible for example that for a d-graph, the set Gp consists of all d- 
dimensional simplices in G. This is still a 0 -graph, a graph with no 
edges because no d-simplex is contained in any other d-simplex. 

Definition 9. Let G be a 2-graph and let / be a locally injective func¬ 
tion on the vertex set of G. Dehne the continuous gradient of / in 
a triangle t = [xyz) with origin x as df = {f{y) — f{x),f{z) — f{x))). 
This assigns to each triangle in G a vector with two real components. 
By taking signs, the vector V/ becomes an element in the hnite vector 
space Z|. This is the discrete gradient on the triangle t rooted at 
the vertex x. 
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Example. 

1) If G is a 2 -graph and x belongs to a triangle t containing two edges 
which both contain x, then Vf{x,t) = (1,1)- Two functions have a 
parallel gradient in t, if and only the sign changes in t happen on the 
same edges of the triangle. 

Remarks; 

1 ) A function / on the vertex set E is a 0 -form. The exterior deriva¬ 
tive df is a function is a 1 -form, a function on the edges of the graph. 
The exterior derivative depends on a choice of orientations on complete 
subgraphs. For 1 forms in particular, where edges have been ordered at 
hrst, the situation at a vertex x dehnes then df{{x,y)) = f{y) — f{x). 
When restricting to a d-simplex, we get d real numbers, we as in the 
continuous forms the gradient. 

2) The ordering of the coefficients of the gradient vector depends on 
the orientation of the triangle. This is similar to the classical case, 
where a triangle in a triangulation of a surface dehnes a normal vector 
at T, once a vertex x of the triangle and an orientation of the triangle 
is given. 

3) The analogue of rank d{f, g) = 2 means that simultaneous sign 
changes happen on one edge of a triangle only, so that the two sign 
change is in the same direction on that edge. Geometrically this im¬ 
plies that the spheres / = 0, 5 ^ = 0 in S'(a;) intersect transversely. Two 
discrete gradients in Z 2 are parallel if and only they are the same be¬ 
cause the only nonzero scalar is 1 . 

The classical Lagrange analysis shows that the critical points of F, the 
place where the Jacobean dF has rank 0 or 1 are candidates for maxima 
of / under the constraint g = c. In the same way as in the continuum, 
we can write down Lagrange equations for any number of functions on 
a d-graph. The most familiar case for two functions: 

Proposition 6 . Let G be a d-graph. If the diserete gradients V/, Vg 
are nowhere parallel, then F = {f = 0, g = 0} is a {d — 2)-graph. 

It should be possible to estimate the measure of the set of global criti¬ 
cal values ii F = (/i,..., fk) are functions on the vertex set of a hnite 
simple graph G. Instead of doing so, we will look later at a Sard setup 
for which the critical values have zero measure. 

Lets look at the example of two functions f,g on a 3-graph. The set 
{F = 0} is the set of simplices, where both /, g change sign. The 
condition dF having maximal rank means that the surface / = 0 on 
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S{x) and the sphere = 0 in S{x) intersect in a union of 1-spheres. 
Here is the Lagrange setup with two functions in three dimensions: 

Proposition 7. If F = (/, g) are two functions on the vertex set of 
a 3-graph and if dF has maximal rank at every x and (c, d) is not in 
F(y), then F = c is either empty of a 1-graph, a finite collection of 
circles. 

Proof. The set {F = c} is a graph whose vertices are the triangles 
(F 3 subgraphs of G) or tetrahedra {K^ subgraphs of G), where both / 
and g change sign. Let a; be a tetrahedron in {F = c}. The maximal 
rank condition prevents parallel gradients like V/ = (1,1,1) = Vg = 
(1,1,1) so that it is impossible to have 3 triangles in {F = c} inside x. 
Assume a single triangle is present where both /, g change sign, then 
there is a common edge, where both /, g change sign and a second 
triangle must also be in {F = c}. We see that there are exactly two 
triangles y, z present on which both /, g change sign. This means the 
vertex x has exactly two neighbors y, z. Each of the triangles y, z has 
two neighbors, the tetrahedra attached to them. We see that {F = c} 
is a graph with the property that every unit sphere is the zero sphere 
Sq. Therefore, it is a hnite collection of circular graphs. □ 

The Lagrange problem extremizing / under the constraint g = cis the 
following: 

Proposition 8 (Lagrange). Given two functions f,g on a 2-graph. 
Extremizing f under the constraint g = c happens on triangles, where 
df and dg are parallel (which includes the case when df = 0 or dg = 0). 
In other words, extrema happen on Lagrange critical points. 

The following maximal rank condition is the same as in the continuum. 
It tells that the graph formed by the simplices having dimension in 
{k,... ,d} on which all functions fk change sign simultaneously is a 
(n — A;)-graph if the rank of the set of gradients in Z 2 is maximal: 

Theorem 9 (Regularity in the commutative setup). Assume G is a 
d-graph and F is a R^-valued function on the vertex set. IfVF has 
maximal rank on every d-simplex, then {F = c} is a {d — k)-graph. 

Proof. For a hxed vertex x and d-simplex X, we have a tangent space 
Zf for which every function fj contributes a vector V/j telling on which 
edges emanating from x inside X, the sign of fi changes. By assump¬ 
tion, the k vectors V/* are linearly independent. 

Use induction with respect to d: if k = d, there is nothing to show 
because by dehnition, the set {F = c} is the set of d-simplices, where 
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all functions change sign and this is a 0-graph as there are no edges. 
We claim that in the unit sphere S{x), the functions fi,..., fk in¬ 
duce k linearly independent gradients V/j on the {d — l)-simplex Y = 
X \x. Indeed, on each triangle, the sum of the dfj values is zero (as 
curl(grad(/)) = 0). A nontrivial relation between the gradients on Y 
would induced a nontrivial additive relation between the gradients on 
X. Now, by induction, the k functions on S{x) dehne a (d — 1 — k)- 
graph. □ 


Remark. One could also try induction with respect to k and consider 
{/i = Cl} which is a (d— l)-graph, by the Sard lemma. The problem is 
that one has to extend / 2 ,..., /^ in such a way on the simplices so that 
one has still maximal rank condition. The problem is that {/i = Ci} 
has now a different vertex set than G and that linking things is difficult. 

Examples: 

1) The case d = 3, k = 2 was discussed in Proposition]^ 

2) In the case d = 4, fc = 2, we want the two functions both to be 
locally injective and the gradients of the two functions /, g not to be 
parallel. The set F = c consists of all tetrahedra K 4 and hypertetra- 
hedra on which both / and g change sign. The gradients restricted 
to the tangent space on S{x) are not parallel and we can apply the 
analysis of the previous case to each unit sphere which shows that in 
S{x) the set F = c is a collection of circular graphs. This shows that 
the unit spheres of F = c have the property that each unit sphere there 
is a circular graph. Therefore F = c is a geometric 2-graph, a surface 
if V/, Vg are nowhere parallel. 

3) In the case d = 4, fc = 3, it is the hrst time that the max¬ 
imal rank condition is not just a parallel condition. Given a ver¬ 
tex X and a tetrahedron t = {x,yi,y 2 ,y 3 )- The discrete gradient is 
(/(di) - f{x),f{y 2 ) - f{x),f{y 3 ) - f ix)). An example of a viola¬ 
tion of the maximal ra nk condition for three functions /, g, h would 
be V/ = (1, 0,1), Xg = (1,1, 0), Xh = (0,1,1). They are pairwise not 
parallel but V/ -|- Xg + Xh = 0. 

4) The non-degeneracy condition is not always needed: for d func¬ 
tions F = ifi,..., fd) on a d-graph, we look at the simplices on which 
all functions change sign. The condition dF 7 ^ 0 implies that the set 
of d-dimensional simplices on which all functions fj change sign are 
isolated. 
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6 . Sard theorem 

Since the set of critical values can hove positive measure if we look 
at the simultaneous solution, we change the setup and look at hyper¬ 
surfaces in hypersurfaces. This will lead to a Sard theorem as in the 
continuum. We will have to pay a prize: the order with which we chose 
the hypersurfaces within hypersurfaces now will matter. 

Definition 10. A function / on the vertex set is called strongly in¬ 
jective if all function values f{xi) are rationally independent. It is 
strongly locally injective if in each complete subgraph, the val¬ 
ues are rationally independent. A list of functions /i,..., /^ is called 
strongly injective if the union of all function values fj{xi) are ratio¬ 
nally independent. 

Strongly injective functions are generic from the measure and Baire 
point of view: given a hnite simple graph G with n vertices. Look at 
the probability space of all functions from the vertex set to [— 1 , 1 ], 
where the probability measure is the product measure on [— 1 , 1 ]'^. 

Lemma 10. For any k, with probability one, a random sample /i,..., /^ 
in is strongly locally injective. 

Proof. There are v = vq + vi + ■ ■ ■ F Vd complete subgraphs in G. 
They dehne the ■ Vi numbers Cij = fi{xj). There is a countably 
many rational independence conditions = 0 to be avoided, 

where Oij are integers. The complement of a countable union of such 
hyperplane sets of zero measure in [— 1 , 1 ]” and consequently has zero 
measure. □ 

Definition 11. Given an ordered list of functions fi,f 2 ,...,fk from 
the vertex list G of a hnite simple graph G = (y,E) and ci,... ,Cfc 
be k values. Let /i = /i. Denote by /2 the function /2 extended to 
{fi = Cl} and by the function /s extended to {fi = ci ,/2 = C 2 } 
etc, always assuming that dehned on = ci,..., fj = cj] does 
not take the value Cj+i. We call the sequence ci,... ,Cfc compatible 
with /i,..., /fc if a sequence f j can be dehned so that none of them are 
not constant. 

Theorem 11 (Discrete Sard for an ordered set). Given k strongly 
injective functions /i,...,/^. For all except a finite set of vectors 
(ci,...,Cfc) the seguence ci,...,Cfc is compatible and the set {fi = 
Cl,..., fk = Ck} is a geometric {d — k)-graph. 

Proof. This follows inductively from the construction. In each step, 
only a hnite set of c values are excluded. □ 
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The assumption is stronger than what we need. It sometimes even 
works in the extreme case of two identical functions: Let / be the func¬ 
tion on the octahedron given the values /i(l) = 13, /i(2) = 15, /i(3) = 
17, /i(4) = 19 on the equator and the value /i(5) = 1 on the north pole 
and the value /i( 6 ) = 31 on the south pole. Lets take ci = 2. Now, 
{/i = Cl} is the cyclic graph with vertices {(51), (512), (52), (523), (53), 
(534), (54), (541)}. The function /2 takes there the values /2(51) = 
(/2(5) + /2(1))/2 = (1 + 13)/2 = 7, 7,(52) = (f,(5) + /2(_2))/2 = 
(1 + 15)/2 = 8, f,{53) = (/2(5) + /2_(3))/2 = (1 + 17)/2 = 9, ^(54) = 
(/2(5) + /2(4))/2 = (l + 19)/2 = 10, /2(512) = (/2(5)+ /2(l)/2(2))/2 = 
(1 + 13 + 15)/3_= 29/3, f,{523) = (^(5) + /2(2)/2(3))/2 = (1 + 15 + 
17)/3 = 35/3, /2(534) = {f,{5) + /2(3)/2(4))/2 = (1 + 17 + 19)/3 = 
37/3, f,{523) = (/ 2 ( 5 ) + / 2 ( 2 )/ 2 ( 3))/2 = (1 + 19 + 13)/3 = 11. Now 
for example, for C 2 = 8.5 the set {/2 = C 2 } is a 0-graph. 

As an example, lets look at the double nodal surface /s = 0 in 
/2 = 0, where /a is the third eigenvector. By Sard, we know: 

Corollary 12. If G is a d-graph and the eigenfunctions / 2 , /s are 
strongly injective not having the value 0, then the double nodal sur¬ 
face is a (d — 2)-surface in G 2 . 

Example: 

1) For the octahedron, the smallest 2-sphere, the spectrum is of the 
Laplacian is 

10,4,4,4,6,6} 

The ground state space is the eigenspace to the eigenvalue 4 and 
three dimensional and spanned by (— 1 , 0 , 0 , 0 , 0 , 1 ), ( 0 , — 1 , 0 , 0 , 1 , 0 ), 
(0, 0, —1,1, 0, 0)}. The eigenspace to the eigenvalue 6 is spanned by 
(1, 0, —1, —1, 0,1), (0,1, —1, —1,1, 0). In both eigenspaces, there are 
injective functions /2 = (—1, —2, —3, 3, 2,1), /s = (1, 2, —3, —3, 2,1) in 
the eigenspace. The graph {/2 = 0} is the cyclic graph C 12 while the 
graph {/a = 0} is CsUCs- The graph {/a = 0} within {/2 = 0} is a two 
point graph. The graph {/2 = 0} within {/a = 0} is not dehned as /2 
extended to the simplex set in a linear way produces a lot of function 
values 0 . 

It leads to a generalization of a result we have shown for 2-graphs: 

Corollary 13. Any compact d-manifold M has a finite triangulation 
which is (d -|- 1)-colorable. 

Proof. By Nash-Tognoli, any compact manifold can be written as a 
variety F = c in some Now just rewrite this in the discrete as the 
zero locus of F = c. □ 
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Corollary 14. The curvature at a vertex x of such a graph triangula¬ 
tion can be written as the expectation of d\ Poincare-Hopf indices. 

If {/ = 0} be the zero locus of f : G ^ R. It would be nice to find a 
smaller homeomorphic graph which represents this set. 

7. Nodal sets 

To illustrate a possible application, lets look at the problem of nodal 
sets of the Laplacian L of a graph. Understanding the Chladni pat¬ 
terns of the Laplacian on a manifold with or without boundary is a 
classical problem in analysis. The nodal region theorem of Courant 
in the discrete also follows from the min-max principle [H |3l |2l] . For 
a general graph G, and an eigenfunction /, one looks at the number 
of connectivity components of = {±/ > 0}. Let Vk be the fc’th 
eigenvector. Since both for compact Riemannian manifolds as well as 
for hnite graphs, the zero eigenvalues are not interesting as harmonic 
functions are locally constant, we look primarily at the second or third 
eigenvalue. The Fiedler nodal theorem assures then that the graph 
generated by Ffc > 0 has maximally k — 1 components. Especially, the 
second eigenvector, the “ground state”, always has exactly two nodal 
components. If an eigenfunction / of the Laplacian is locally injective 
and has no roots, we can look at its nodal surfaces / = 0 separating the 
nodal regions. I / does not take the value 0, then {/ = 0} is defined 
even if / is not locally injective. As in the continuum, one can ask how 
big the set {n G U | f{v) = 0} can become. 

If G is a 2-sphere, then the two nodal surface is a simple closed curve 
in G. How common is the situation that the ground state does not 
take the value 0 and is locally injective? The situation that 0 is in the 
range of the ground state appears to be rare. For random 2-spheres (of 
the order of 500 vertices generated by random edge rehnements from 
platonic and Archimedean solids) we get a typical ground state energy 
in the order of 0.08 and the third eigenvalue in the order 0 . 2 . 

[25] note the following eigenfunction principle: any eigenfunction 
to an eigenvalue 0 < A < n takes the value 0 on every vertex of 
degree n — 1, if n is the number of vertices: Proof: let / be the 
eigenfunction to an eigenvalue A. The function / is perpendicular 
to the harmonic constant function so that YIxgv /(^) “ From 
Lf{v) = {n- l)f{v) - we get nf{v) = Xf{v) 

which by assumption implies f{v) = 0 . 
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Example. 

For a wheel graph with n vertices, there is one eigenvalue n with eigen¬ 
vector (1 — n, 1,1 ,..., 1), an eigenvalue 0 with eigenvector (1,1 ,..., 1). 

All eigenfunctions to eigenvalues between take the zero value some¬ 
where. For example, in the case n = 7, the eigenvalues are {0, 2, 2,4,4,5, 7} 
with eigenvectors [0, —1, —1, 0,1,1, 0], [0,1, 0, —1, —1, 0,1], [0, —1,1, 0, —1,1, 
[0,-1, 0,1,-1,0,1], [0,-1,1,-1,1,-1,1] and [-6,1,1,1,1,1,1], 

Let G be a d-graph. Let / be the ground state, the eigenvector to the 
hrst nonzero eigenvalue A, the spectral gap. Let d denote the exterior 
derivative. There will be no confusion with d also denoting the dimen¬ 
sion of G. We assume that the eigenvalue A is simple, that / has no 
roots and that df has no roots. This assures that Z = 0 is a geometric 
(d — l)-graph. As L = dd* on 0 forms, we get from d*df = \f that 
dd*df = Xdf so that df is an eigenfunction to the one form Laplacian 
Li = dd* + d*d. The set {df = 0} of all triangles and simplices where 
both / and df changes sign is the same than / = 0. 

What is the topology of the hypersurface Z 2 = {v 2 = 0}? For a 2- 
sphere G, we know that there are two components so that the principal 
nodal curve Z 2 has to be a circle. Is the nodal curve to the second 
eigenvalue a d-sphere, if G is a d-graph? We believe that the answer 
is yes and robust. If / should have roots, we can add a small random 
function with \g{x)\ < e. We expect that for sufficiently small e > 0 
and almost all g, the surface Z = [f + g = 0 } has the same topology. 

Does the topology of the nodal manifold to the second eigenvalue de¬ 
pend on the topology of G only? 

We think the answer is yes as Z has to be a connected surface and 
that going from a genus k to a genus k + 1 surface can not happen so 
easily. More generally, we expect: if G and H are homotopic graphs of 
the same dimension, then the second nodal manifolds of G and H are 
homotopic. 



22 


OLIVER KNILL 



Figure 3. Level surfaces of the second eigenfunction 
/2 of the Laplacian on G 2 , where G is a 2-sphere. The 
dividing surface is then always a Jordan curve. 
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Figure 4. Level surfaces of the second eigenfunction 
of the Laplacian, where G are 3-spheres. The dividing 
surface is then a 2-graph. We believe that {/2 = 0} 
is always a 2-sphere, if G is a 3-sphere and /2 is the 
eigenvector to the smallest nonzero eigenvalue. A more 
general question is whether two whether the nodal region 
{/2 > 0} is always simply connected if G is a d-sphere. 
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Figure 5. Level surfaces of eigenfunctions of the 
Laplacian on a graph Gi, where G is the icosahedron. 
The Courant-Fiedler nodal theorem tells that the num¬ 
ber of positive nodal regions of is < n. It is conhrmed 
in the pictures seen. The first picture shows the second 
eigenvector. The first one is constant and has only one 
region. 




A SARD THEOREM FOR GRAPH THEORY 


25 





Figure 6 . Level surfaces of eigenfunctions of the 
Laplacian on G 2 , where G is the icosahedron. Also this 
illustrates the Courant-Fielder theorem. The hrst pic¬ 
ture shows the second eigenvector for which the surface 
Z 2 divides the graph in to two regions. The hrst one is 
constant and has only one region. 
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8. Summary 

A d-graph is a finite simple graph for which every unit sphere is a {d—1) 
sphere. Given k locally injective real-valued functions fi,..., fk on the 
vertex set of a d-graph, we can dehne the zero locus Z = {F = 0} 
of F in G as the graph (V, E) with vertex set V consisting of com¬ 
plete subgraphs in G, on which all the functions fj change sign and 
for which two vertices are connected if one is a subgraph of the other. 
The graph Z is a subgraph of the barycentric rehnement Gi but like 
varieties in the continuum, the graph might be singular and not be 
a (d — A;)-graph. We can dehne a discrete tangent bundle on Gi so 
that the projection of V/ in T^G = to T^H is zero, establishing 
the graph analogue of “gradients are perpendicular to level surfaces”. 
If every unit sphere S{x) in Gi has the property that Z H S{x) is a 
{d — k — l)-sphere in the (d — l)-sphere S{x), then F = c is a (d — /c)- 
graph H. The later plays the role of {/i = Ci, • • • = fk = Ck} for k 
differentiable functions / on a smooth manifold M for which the max¬ 
imal rank condition rank(dF)(a;) = d — /c for F = {fi,..., fk),x G Z 
is satished. Discrete Lagrange equations help to maximize or minimize 
/ under a constraint g = c: as in the continuum, critical points need 
parallel gradients Vf^Vg. 

If the ordered zero locus is dehned by looking at hypersurfaces on 
successive barycentric rehnements, we never run into singularities if 
the functions are strongly locally injective in the sense that the func¬ 
tion values of fj on the vertex set are rationally independent on each 
complete subgraph. The Sard theorem states then that the graph 
{/i = Cl,... ,/fc = Ck} is a (d — A;)-graph inside the d-graph Gk for 
all c not in a finite set. An application is the observation that any d- 
dimensional projective algebraic set admits an approximation by mani¬ 
folds with a triangulation which is a minimally {d+ l)-colorable graph, 
using an explicitly constructed d-graph determined by the equations 
dehning V. 

The possibility to dehne level surfaces in discrete setups can be illus¬ 
trated in the case of eigenfunctions of the Laplacian. Pioneered by 
Chladni, the geometry of nodal surfaces is important for the physics 
of networks. Of special interest is the ground state /2 of a com¬ 
pact Riemannian manifolds or hnite graph. The nodal surface /2 = 0 
is a (d — l)-dimensional hypersurface and the double nodal surface 
/2 = 0, /a = 0 is generically a codimension 2 manifold. If the manifold 
or graph is a 3-sphere, the nodal surface is two dimensional and the 
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double nodal surface is a collection of closed curves in S^. We asked 
here whether /2 = 0 have positive genus and whether /2 = 0, /s = 0 
can be knotted. Graph theory allows to investigate this experimentally. 
The answers for Riemannian manifolds and graphs are expected to be 
the same. 

Lets look at some multi-variable calculus termiuology in a 2-graph 
translated to graph theory: 


Critical point x 

Sj {x) is not contractible 

Discriminant D 

PoincarCHopf index if{x) = 1 — x{Sj{x)) 

D>0Ja:,<0 

if{x) = 1, Sj{x) = Sf{x) 

> 0, > 0 

if{x) = 1, Sj{x) = 0 

D <0 

if {x) < 0 

D = 0 

not locally injective 

Level curve f{x, y) = 0 

zero locus / = 0 in Gi. 

T.M 

maximal simplex t in G containing x 

Tangent vector T^M 

V/ = sign(/(i/) - f{x)) \ yet} 

Lagrange equations 

V/ = XWg or Vg = (0,0) 



Figure 7. The nodal hypersurface {/ = 0} of 
the ground state. The length of the surface is 20, 
the two nodal domains have 30 triangles each. The 
Cheeger number (see [1] for a discrete version), c(/) = 
|C'(/)|/min(|yl(/)||R(/)|) satsihes the Cheeger inequal¬ 
ity. 
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Figure 8. The spectrum of refinements G,Gi,G 2 ,G 3 
of the icosahedron graph G. We noticed that the spec¬ 
tral integrated density of states converges to a function 
which only depends on the dimension of the maximal 
complete subgraph. While for graphs without triangles, 
the convergence of Fg{x) = X\j^x\ to the limiting func¬ 
tion 4sin^(7ra;) is immediate, in higher dimensions, the 
limiting function appears to be non smooth. See HZ]. 
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9. Questions 

A) In the commutative case, where F = c is a subgraph of the barycen- 
tric rehnement Gi, the set of critical values can have positive measure. 
It would be nice to hud an upper bound on the measure of constants c 
for which the graph {F = c} C Gi can be singular. 

B) In the context of graph colorings, the following question is analog 
to the Nash embedding problem: Is every d-graph a subgraph of 
a barycentric rehnement G„ of some geometric d-graph? If the answer 
were yes, we would have a bound d + 1 for the chromatic number of G. 
A special case d = 3 would prove the 4-color theorem. We approached 
this problem by writing the sphere as embedded in a 3 sphere then 
making homotopy deformations to render the sphere Eulerian and so 
4 colorable coloring in turn the embedded sphere. While a Whitney 
embedding of a graph is possible if we aim for a homeomorphic image, 
we don’t have yet a sharp discrete analogue of the classical Whitney 
embedding theorem, hxing the dimension. Realizing a graph as a sub¬ 
graph of a product of linear graphs is analogue to an isometric embed¬ 
ding of a Riemannian manifold in some Euclidean space is a discrete 
Nash problem. 

C) If f,g are two strongly locally injective function on some d-graph. 
Under which conditions is it true that the level set g = d in f = c 
is topologically equivalent to the level set g = d in f = c7 Random 
examples in a 3-sphere show that the answer is no in general. It might 
therefore be possible that /2 = 0 in /a = 0 is different that /a = 0 in 
/2 = 0 but we have not yet an example, where this difference takes 
place. 

D) Assume G is a 3-sphere and the ground state fi (eigenvector to the 
smallest positive eigenvalue) does not take the value 0. Is the nodal 
hyper surface {/ = 0} always a 2-sphere? The case where 0 is in the 
image of the ground state fi only appears in very rare cases. As long as 
the O’s are isolated, we can randomly change the value on such places 
and not change the topology of the nodal surface {/ = 0}. 
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